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Abstract 



The non-abelian generalization of the Born-Infeld non-linear lagrangian is ex- 
tended to the non-commutative geometry of matrices on a manifold. In this case 
not only the usual SU(n) gauge fields appear, but also a natural generalization of 
the multiplet of scalar Higgs fields, with the double-well potential as a first approx- 
imation. 

The matrix realization of non-commutative geometry provides a natural frame- 
work in which the notion of a determinant can be easily generalized and used as the 
lowest-order term in a gravitational lagrangian of a new kind. As a result, we obtain 
a Born-Infeld-like lagrangian as a root of sufficiently high order of a combination of 
metric, gauge potentials and the scalar field interactions. 

We then analyze the behavior of cosmological models based on this lagrangian. 
It leads to primordial inflation with varying speed, with possibility of early decel- 
eration ruled by the relative strength of the Higgs field. 
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1 Introduction 



1.1 The origin of the Born-Infeld lagrangian 

Among numerous cosmological models using scalar field as a source of primordial cosmic 
energy and the subsequent inflation there is usually no limit on the field strength: as 
a matter of fact, the scalar field $ can take on arbitrarily high values. This is true 
in particular for the simplest model with V(0) = X(j) 2 (<p 2 — 7 2 ) "double well" potential. 
However, one may point out that principles of quantum field theory, when combined with 
General Relativity, should lead to the existence of a cut-off for any field's value, because 
at certain intensity the very notion of space and time will not be valid anymore. More 
precisely, it is because Heisenberg's uncertainty principle induces a spontaneous birth of 
black holes, when applied to the gravitational field. 

This is why it seems reasonable to investigate theories in which such a cut-off is 
incorporated from the very beginning. First such attempt concerned exclusively the 
electric field, whose infinite value at r — > in the Coulomb law should be avoided in 
order to keep the energy finite. 

To this purpose, G. Mie introduced the notion of maximal field strength, E , in 
order to make it impossible for any electric field to go beyond this value. He modified 
Maxwell's theory by introducing the following non-linear lagrangian density for pure 
electric field: 



E 2 



Although the non-linear theory derived from this lagrangian enabled Mie to obtain a 
non-singular, finite energy solution, it was clear that such a lagrangian can not represent 
a Lorentz-invariant theory, especially that from the beginning there was no magnetic field 
contribution at all. This is why Born and Infeld [2] introduced a Lorentz-invariant la- 
grangian density (when constructed on a pseudo-Riemannian manifold, it is also invariant 
under diffeomorphisms), defined as follows: 

Cbi(9,F) = 




(2) 

The constant (3 appears for dimensional reasons and plays the role of the upper limit of 
the electric field in Mie's non-linear electrodynamics. Defining 

P=\f^ v F» v and S = ^ F^ , with F^ = ± e» vXp F Xp 

the Born-Infeld lagrangian can be recast in a simple form: 

Cbi = (3 2 1-VI + 2P-S 2 (3) 
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The idea of a non-abelian generalization of Born-Infeld theory lagrangian has been in 
the air already at the end of the seventies. Hagiwara has discussed various possibilities 
in however, he did not try to find soliton-like solutions. In 1997 Tseytlin f I argued in 
favour of the symmetrized trace prescription which reproduced in the first 4 orders the 
string effective action for gauge potential. Finally Park introduced yet another non- 
abelian generalization and investigated qualitative behavior of instanton-like solutions. 

It is interesting to note that the non-commutative generalizations of geometry of the 
space-time lead quite naturally to this class of lagrangians. 

1.2 Arguments in favour of non-commutative geometry 

Simultaneous consideration of the two most important new physical theories of this cen- 
tury, General Relativity and Quantum Mechanics, have not yet produced a common tool 
for the description of the nature of space-time at the microscopic level. The General Rel- 
ativity develops our knowledge about global properties of space and time at very large 
distances, and raises questions concerning the global topology of the Universe. 

The methods of Differential Geometry, the best adapted mathematical language of 
this theory, are very different from the methods of Quantum Physics, in which one stud- 
ies the properties of the algebra of observables, considering the state vectors, as well as 
geometric points and trajectories, as artefacts and secondary notions. This approach has 
been inspired by the works of John von Neumann [9], and has much in common with the 
non-commutative geometry, where the very notion of a point loses its meaning. 

There are several well-known arguments in favour of the point of view according to 
which the dynamical interplay between Quantum Theory and Gravitation should lead to 
a non-commutative version of space-time. Let us recall the most frequently cited ones: 

* A semi-classical argument that involves the black hole creation at very small 
distances: as a matter of fact, if General Relativity remains valid at the Planck scale, 
then any localization of events should become impossible at the distances of the order 

of Xp = Indeed, according to quantum mechanical principles, to localize an event 

in space-time within the radius A / ~ a, one needs to employ energies of the order of 
a -1 . When a becomes very small, the creation of mini black holes becomes possible, 
thus excluding from the observation that portion of the space-time and making further 
localization meaningless. Therefore, the localization is possible only if we impose the 
following limitation on the time interval: 

Ax°(SAx fc ) > X 2 P and Ax k Ax m >X 2 P . (4) 

in order to avoid the black hole creation at the microscopic level. 

** The topology of the space-time should be sensitive to the states of the fields 
which are in presence - and vice versa, quantum evolution of any field, including gravity, 
should take into account all possible field configurations, also corresponding to the fields 
existing in space-times with radically different topologies (a creation of a black hole is 
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but the simplest example; one should also take into account other "exotic" configurations, 
such as multiple Einstein- Rosen bridges (the so-called " wormholes") , leading in the limit 
of great N to the space-time foam). 

Now, as any quantum measurement process may also lead to topological modifica- 
tions, again the coordinates of an event found before and after any measurement can not 
be compared anymore because they might refer to uncompatible coordinate patches in 
different local maps. As a result, quantum measures of coordinates themselves become 
non-commutative, and the algebra of functions on the space-time, supposed to contain 
also all possible local coordinates, must be replaced by its non-commutative extension, 
better adapted to describe the space-time foam. 

*** Since the coordinates are endowed with a length scale, the metric must enter 
at certain stage in order to measure it. After quantization, the components of the tensor 
g^y become a set of dynamical fields, whose behavior is determined by the propagators 
and, at least at the lowest perturbative level, by two-point correlation functions. As any 
other field, the components of the metric tensor will display quantum fluctuations, making 
impossible precise measurements of distances, and therefore, any precise definition of 
coordinates. Let us now consider possible consequences of the above assumptions. 

2 The Non- commutative Born-Infeld lagrangian 

2.1 Non-Commutative Einstein-Hilbert lagrangian 

Let us assume therefore that at the Planck scale not only the positions and momenta 
do not commute anymore, but also the coordinates themselves should belong to a non- 
commutative algebra. As a direct consequence, all functions of coordinates that served 
to describe various geometrical quantities should also belong to this algebra, although 
some of them can belong to its center and commute between themselves. 

The non-commutativity of coordinates means also that the notion of distance between 
two points of space-time can not remain symmetric: if the quantum indetermination 
principle applies to the measurements of positions, then there is no guarantee that the 
distance mesured first from a point A to its neighbor B will remain strictly the same 
when we try to measure next the distance from B to A. This amounts to the introduction 
of a non-symmetric metric tensor, g^ ^ g^. In his efforts to create a general field 
theory unifying gravity with other fields Einstein considered such a possibility in his 
late years, but without quantum theoretical motivation, and his investigations remained 
purely classical. 

Now, combining two direct consequences of non-commutative character of space-time 
coordinates and their functions leads quite naturally to the following general form of the 
metric tensor: 

9v = 9i i + Gi S + f% +4, (5) 
In this formula we have separated the commutative and non-commutative parts, as well 
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as the symmetric and antisymmetric parts: the terms g° and are supposed to belong 
to the center of our hypothetic non-commutative algebra replacing the usual algebra of 
smooth functions of coordinates, with g° = and = — /&. The terms denoted 
by capital letters correspond to the non-commutative operators with similar symmetry 
properties of their tensor indices, 

Gij = Gji and Fy = — (6) 

Assuming that both non-standard effects, the asymmetry of the metric tensor and the 
non-commutative character of its components are coming from the same quantum source 
whose manifestations become detectable only on Planck's scale, it is natural to assume 
that the commutative part of is symmetric, whether the skew-symmetric part is also 
a genuine non-commutative operator. This leads us to the conclusion that it is enough 
to keep only two terms in the expression [HJ 

9ij=9y + fij, (7) 

Here g° also belongs to the non-commutative algebra of operators, but it remains in its 
center and should commute with everything else. Having this point made clear, we shall 
drop the distinctive upper signs and use the simplified notation + Fy, remembering 
the particular character of each term. 

Thus the quantum theoretical arguments lead us naturally to the generalized form 
of metric that served as the basic ingredient in the Born-Infeld generalization of electro- 
magnetic lagrangian. 

The generalized metric tensor should be now used in order to produce a new vari- 
ational principle defining the dynamical behavior of fields. We should follow the usual 
scheme which consists in constructing the series of invariants of metric tensor starting 
from the lowest order which is the volume element. In classical differential geometry one 
has: 



5S = 5 



(8) 



where g denotes the determinant of g^ so that yfg defines the volume element in local 
basis; R is the Riemann scalar, the next term is the Gauss-Bonnet invariant, and so forth. 
The coefficients A, 1/G and 7 define the relative strength of each contribution. We see 
now that the first term in the integrand has the form of the Born-Infeld lagrangian, 
and is probably a good candidate for Planckian generalization of Einstein's cosmological 
term. Taking into account the relative strengths of known classical fields, it is clear that 
the skew-symmetric tensor Fy appearing in the generalized metric does not describe the 
classical electromagnetic field, but rather the quantum field effects responsible for the 
early behavior of the Universe right after the Big-Bang. 

As in the Kaluza-Klein type theories, one can also expect the contribution of the 
F-field coming from the next term containing the Riemann scalar; however, they enter 
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only through their first and second derivatives, as all other componenets of the metric 
tensor, therefore they will alter only the derivative part of the lagrangian, and not the 
scalar potential part found in the cosmological term. 

In what follows, we shall expose in a concise way, on the example of the simplest finite 
non-commutative algebra, which is the algebra of complex nxn matrices, how almost all 
the notions of usual differential geometry can be extended to the non-commutative case. 
We shall also show how the gauge theories and the analogs of the fibre bundle spaces and 
Kaluza-Klein geometries can be generalized in the non-commutative setting. The finite 
version of non-commutative geometry is also the best adapted for the generalization of the 
determinant of metric tensor needed for the construction of the Born-Infeld lagrangian. 

2.2 The Non- Commutative Matrix Geometry 

We shall generalize now the non-commutative Maxwell theory developed in [16J in order 
to obtain a Born-Infeld like theory. Let us resume the notations and language of the 
theory. We consider the algebra A = C°°(V) <g> M n ( C) with the "vector fields" spanned 
by the derivations of C°°(V) and inner derivations of M n (C). The differential algebra is 
generated by the basis of linear 1-forms acting on the derivations. In order to construct 
a gauge theory, we must consider a finite projective module over A. By analogy with the 
Maxwell theory, we will consider the simplest one, i.e. the free module of rank 1 over 
A, which can be identified with A itself. Then one defines a gauge by the choice of a 
unitary element e of A, satisfying h(e, e) = 1, with h an hermitian structure on A. Then 
any element of A can be written in the form em with m G A and a connection on A is 
a map: 

V : .4 ^(.4), e m i-> (Ve) m + e dm (9) 
In the gauge e, the connection can be completly characterized by an element lj oiD, 1 (A): 

Ve — e to . 

One can also decompose lj in vertical and horizontal parts: 

lj = uJh + lj v with LJh = A, lj v = 6 + (p (10) 

Here A is like the Yang-Mills connection, whereas 9 is the canonical 1-form of the matrix 
algebra, and plays the role of a preferred origin in the affine space of vertical connections. 
It satisfies the equation: 

do + e 2 = o 

Then is a tensorial form and can be identified with scalar field multiplet. 

Choosing a local basis of derivations of A: {e M , e a }, where for convenience e^ are outer 
derivations of C°°(V), and e a = ad(X a ), with {A a } a basis of anti-hermitian matrices of 



5 



M n (C), are inner derivations. 

The dual basis will be denoted by {9^,9 a }. In this particular basis, we have: 

a = , e = -x a e a , <p = cj) a e a 

If we choose the connection to be anti-hermitian, we can write <f) = <j) b a \fy8 a . The curvature 
tensor associated with uj is : 

Vl = duo + uo 2 

we can also define the field strength: 

F = dA + A 2 . 

Then by "dimensional reduction" one can identify: 

Qa^ = -D^a Slab = [4>a, 4>b] ~ C^c 

where C^ b are the constant structure in the {A a } basis. 

A gauge transformation is performed by the choice of a unitary element U of M n (C), 
satisfying h(eU, elf) = 1. Then in the gauge e' = eU 

J = U~ x uU + U~ 1 dU 

6 is invariant under gauge transformations, then 

A 1 = U- l AU + U^dU, <f)' = U'^U 

2.3 The non-commutative version of the Born-Infeld lagrangian 

In this section, we will essentially recall the non-commutative generalization of the Born- 
Infeld lagrangian proposed in By analogy with the abelian case, we want the la- 
grangian to satisfy the following properties: 

1) One should find the usual NC Yang-Mills-Higgs theory in the limit j3 — ► oo 

2) The (non-abelian) analogue of the "electric" field strength should be bounded from 
above when the "magnetic" components vanish. 

(To satisfy this particular constraint, we must ensure that the polynomial expression 
under the root sign should start with terms 1 — (3~ 2 {E a ) 2 + ... when B a = ) 

3) The action should be invariant under the action of the automorphisms of A (it 
includes diffeomorphisms and gauge transformations). 

4) The action must be real. 
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This enables us to introduce the following generalization of the Born-Infeld lagrangian 
density for a non-commutative gauge field: 

y/det\g\ -{|det(l®# + J®(l\} 1/in (11) 

and Q = VL a pL al3 with L a P the generators of the fundamental representation of SO (A + 
n 2 — 1), Q a p are the components of the curvature defined in previous section, and then 
are anti-hermitian elements of M(n, C). In the expression above, J denotes a unitary 
2x2 matrix satisfying J 2 = — 12, thus introducing a quasi-complex structure. This extra 
doubling of tensor space is necessary in order to ensure that the resulting lagrangian is 
real. At the same time, it raises the degree of the polynomial under the root up to An. 
We are left with the root of order An, so that the invariance of our action under the 
space-time diffeomorphism is preserve. 

Let us recall a few arguments in favour of this construction: 

The simplest way to generalize the Born-Infeld action principle in the framework 
of the non-commutative matrix geometry seems at first glance the substitution of real 
numbers by corresponding hermitian operators, like in quantum mechanics. Then one 
would arrive at the following expression: 

where l n and iT a are n x n hermitian matrices. What remains now to make the gen- 
eralization complete, is to extend the notion of determinant taken over the space-time 
indeces in the usual case, i.e. the determinant of a 4 x 4 matrix, to a notion of deter- 
minant taken in the tensor product of derivations and matrix indeces of the algebra A. 
Then one would replace the objects in (pj) following the procedures in (fl"2|) but it leads 
to a complex lagrangian. In order to avoid this problem, it is necessary to start from an 
other form of the usual Born-Infeld lagrangian: 



Sbi[F, 9] = (3 2 (j\g\ - \det c (l 2 ® 9^ + /T 1 J ® iF^ ) | % \ d A x , (13) 



where J is a 2 x 2 unitary matrix whose square is equal to — 12. With the correspondence 
displayed in (|T2"1) . we arrive at the action principle displayed in (fTT| which satisfies all 
the requirements we asked for, 1), 2), 3) and 4). 



The lagrangian in (fTTJ) contains the contribution of two types of fields: the classical 
Yang-Mills potential, A = A^O 1 *, corresponding to the usual space-time components of 
the connection one-form, and the scalar multiplet coming from its matrix components 
(ft = (p a 6 a = (f) b a \b9 a . In the case when = 0, this lagrangian coincide with the one studied 
in [Zj. For cosmological considerations, we willrestrict ourselves to a qualitative analysis 
of the case when the space time components of Q do vanish F^ v = 0, leaving only the 
contribution of scalar multiplet degrees of freedom. 
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2.4 The reduced lagrangian for scalar fields 



Let us recall the notations which will be used in the subsequent calculations. The basis 
of matrix representation of the SU (2)-algebra is chosen as follows: 



A a — — i(Ja A a A;, — — 5 a b + e a b c A c [A a , A&] — C c ah — 2e a b c \ c 



(14) 



Now we have to evaluate the determinant of the following matrix: 



1 iD<j> 
-iD6 1 + iH 



where 

I J a,6=l,2,3 ^ 

From now on, we choose the simplest ansatz 1 with one scalar field only: 

= if e 

After some algebra, we get the following result: 



(15) 



(16) 



L = 1 - {1 + 6(D<p) 2 + 9(LV) 4 + I6p 2 (<p - I) 2 } 1 v/1 + V(y>- 1) 



(17) 



In [Zj we have introduced a new non-abelian generalization of the Born-Infeld la- 
grangian, and found a family of non-singular soliton-like solutions, using 't Hooft's ansatz 
for the SU{2) gauge potential. In [6j, we have generalized this approach in order to in- 
clude scalar multiplets arising naturally in the non-commutative geometry of matrix 
valued functions. In the case when all degrees of freedom are reduced to a single scalar 
field cp, we have investigated homogeneous time dependent solutions. As expected, the 
regular trajectories in the phase space ((p, ip) turn out to be comprised in a finite domain, 
thus confirming the existence of a finite bound on field strength. This particular behavior 
of the Born-Infeld like scalar field makes us believe that it can be used in the framework 
of cosmological models on a footing similar to the so called inflaton field. 

This lagrangian will become the first ingredient for the cosmological model we want 
to introduce. The next termin the lagrangian will be the usual Riemann scalar term 
providing the minimal interaction with the scalar Born-Infeld field through covariant 
derivatives and through the direct coupling with the symmetric part of metric tensor and 
interactions with the Yang-Mills fields. 

X A detail analysis of other possible ansatzs is performed in [H] 
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3 Non-commutrative Born-Infeld cosmology 



3.1 General considerations 

The study of Born-Infeld lagrangians both in the abelian [Hj and the non-abelian [TH] 
case or as condensates [12J has been widely developed in the last few years in relation 
with string theory, where Born-Infeld type lagrangians appear in a natural way [T3| IT1]. 
Incidentally, scalar fields can be introduced in diffent ways in this context by considering 
the so-called Dirac-Born-Infeld (DBI) action wich could be relate to the usual Born-Infeld 
action by T-duality and dimensional reduction. Then the action we are considering can 
be seen as a non-commutative matrix DBI action. 

It has been argued that the U(l) case shows the problem of the axial symmetry 
induced by the electromagnetic field F^ u which forbids homogeneous space-time solutions. 
On the other side in the framework of a non-abelian gauge theory spherical symmetry 
can be restored and a cosmological model in the standard FRW metric can be analyzed 
[m ESI- m this case some interesting results have been obtained. In fact it has been 
demonstrated that this kind of lagrangian can furnish an inflationary scheme [10] or 
provide a theoretical framework to the phantom field models ^2] introduced in several 
theoretical attempts to explain the dark energy [T8l fTO} 120 } l2T]. 

The introduction of a scalar multiplet, reduced to one |6j, inside the SU{2) gauge group, 
goes in this direction, so that in the following we will consider a standard FRW metric 
as the gravitational framework of our cosmological model. One may think that this kind 
of scalar field could drive the first stages of the evolution of the universe. 

The most widely accepted scheme is to consider a scalar field rolling down in its po- 
tential well slowly enough to provide a cosmological constant like term in a first phase, 
driving an exponential expansion. In a second phase the scalar field decays in a true 
vacuum configuration oscillating coherently around the minimum of its potential and 
transferring its energy to the standard matter fields (reheating). Many approaches have 
been proposed in time (old inflation, new inflation, (2jj, chaotic inflation [28j, see [201 for 
a review) also with completely different mechanism (Starobinsky scalaron |30j). 
The Born-Infeld scalar field may be interpreted just as an effective bosonic field conden- 
sate and not a fundamental field. 

In our model we shall consider minimal coupling with gravity, obtained by adding the 
usual Einstein-Hibert lagrangian of gravitational field, and by replacing in Euler-Lagrange 
equations all derivatives by their covariant counterparts. The lagrangian density will be 

£ = £n C &i(0, <9 M 0, fiv) + ^q R ■ ( 18 ) 

The Born-Infeld type scalar field lagrangian L nc u contains two parameters: the maximal 
value of the field strength (3 and the characteristic mass 7 of the scalar field. Together 
with Newton's constant G, our problem will contain 3 explicit parameters with dimension 
of mass. By letting G — > the gravitational interaction is decoupled and we recover the 
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model of scalar Born-Infeld field studied recently in jH]. We can also recover the usual 
</> 4 theory when the parameter f3 tends to infinity. In what follows, a systematic study of 
the dependence of the cosmological model on the three aforementioned parameters will 
be presented. 

3.2 Cosmological model with Born-Infeld scalar field 

The lagrangian fTTjl can be write after a rescaling of the scalar field: 

1/4 



■"ncbi 




(l-r 2 2 (t)) 2 + y/T 2 (-V37 + 0W 



(19) 



where the parameter (3 is the analog of the BI cutoff for the non linear electromagnetic 
theory while 7 gives account of the mass-shell of the scalar field and depends from the 
non commutative algebra considered for the gauge group. It is obvious that such a model 
contains in itself the 4 theory 2 so we expect to obtain a standard behavior when the 
Born-Infeld parameter (3 tends to infinity. To implement this scheme in presence of 
gravity, the minimal coupling with the Hilbert-Einstein lagrangian can be considered, 
with the action given in the Eq. (fT8|) . It reads explicitly as follows: 



The two functions T [</>(£), </>(£)] and II[0(£)] are defined as follows: 

X [</>(£),</>(£)] = ((l-r 2 2 W) 2 + yn 2 [(0(£)]) (21) 

I1[0(£)] = /rV(i)(0(i)-V3 7 ). (22) 
Varying (|271j) with respect to the metric we get two cosmological equations: 



" Ht) + -L = \K- |r 2 2 (t)(i-r 2 2 (t))T-i[0(£),0(£)] fi + ^n 2 [0(£)]V /2 + 



a 2 (£) a 2 (£) 3m 2 Pj f ^ wv ^ ^ v " Ln/,rwj \ 9 



4 



1/2 



T4[0(£),0(£)](l + -n 2 [0(£)]) -\) , (23) 



2 this can be observed while developing the expression iflflll in a Taylor series as it will be shown in 
the next section 
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which is the energy equation (Friedmannn equation) and 

which is the pressure equation. 

The third equation is obtained varying the lagrangian with respect to the scalar field. 
It can be shown that it coincides with the conservation equation for the energy momentum 
tensor, namely the Bianchi identity [HI]. We have: 

4>(t) (1 + V[0(t)]) | (1 - /r 2 2 (*)) 2 + yn 2 [0(t)] (1 - 3/r 2 2 (t)) } + 

- v / 3 7 )(20(t) - \/3 7 ) (T[0(t),0(t)]) x (25) 

x (3 - /rV 2 (t) + V[0(t)]) - 6/?- 2 2 (t) (1 - /rV(*)) (1 + ^{m) + 



+3^|| m (l + ^im) (l - /5~ 2 2 (*)) (l " 2/?" 2 2 (t) + ^*(t) + yH 2 [0(t)]) . 

From the lagrangian or directly from the cosmological Eq. (J2Bj) and Eq.lj2%|) we find that 
the Higgs-type scalar field is characterized by its energy density and the pressure term: 



Pncu = P 2 j/r 2 2 (*)(i - /3- 2 2 (t))T-t[0(t),0(t)] (1 + ±iPW)]f + ^W)A{t)\ - 1 } , 

(26) 

Pn C M = /? 2 {l-Ti[0(t),0(t)]} , (27) 
which allows us to write the barotropic factor as well 

Pncbi /oo\ 

Wncbi= • (28) 

Pncbi 

At the first sight the system appears quite complex. This is why, as a first step in 
what follows, we will perform a qualitative analysis of the model through numerical study 
of the solutions. In this way it will be possible to get significant conclusions about the 
dynamics of the system and its cosmological interpretation. 
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3.3 The Born-Infeld relaxed limit 

We have seen that in the general case the NCBI lagrangian (TT91) is characterized by 
means of two parameters like in the standard Higgs theory. However in this case their 
origin is strictly connected with the non-abelian character of the gauge field and with the 
Born-Infeld model itself. 

One of these parameters, namely 7, comes from the non abelian character of the Yang- 
Mills theory and it is connected with the scaling of the non commutative algebra. 0[E|. 
On the other side, the parameter (3 arises in relation to the choice of a Born-Infeld type 
lagrangian for the gauge theory [H]. It acts as the Born-Infeld typical threshold for the 
field intensity. The two parameters have different dimensions. In our formalism the NC 
Yang-Mills gauge field fl^ which also contains the scalar field El has the dimensions 
of [cm] -1 as the scalar field itself and 7. This result indicates that one can consider 7 
as a mass-shell parameter of the scalar field imitating the typical mass parameter of the 
Higgs field. Regarding (3 it is easy to check that (TH^I it has to be [cm] -2 . 

Now, to compare with the linear limit of the NCBI lagrangian, we can develop it in 
series of (3. As expected, the NCBI model shows, in the linear limit, a typical Higgs 
behavior with a "top hat" potential: 



The shape of the potential depends only on the coefficient 7, which sets of the mass of 
the scalar field. Of course the analogy with the Higgs field is not complete because of 
the Born-Infeld linear limit has no coupling constant to compare with the usual Higgs 
sector. The different choices of the values of parameters f3 and 7 will generate not only 
a whole spetrum of inflationary behaviors of our model, but will also create a possibility 
for the acceleration as well. 

4 The phase space study 
4.1 Preliminary Considerations 

The field equations of the NCBI model appear extremely complicated and there is little 
hope to solve them analytically. But qualitative analysis can be performed instead. To do 
this we should analyze the phase space of the system considering simultaneous evolution 
of the scale factor and of the scalar field. We assume a flat space-time (k = in Eq. 
p^jl -pl jl ). since more and more results coming from the CMBR observations seem to 
confirm the spatially flat model. It also seems that this assumption is well supported by 
other astrophysical data. After the COBE results of 1992 [23] other new data have been 
obtained in the last years. The balloon-based experiments [2E] and the 2003 the WMAP 
satellite provided stringent estimates of the last scattering surface anisotropies indicating 
that the actual spatial curvature is very close to zero. 




(29) 
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A few more considerations will enable us to better understand the physical properties 
of the NCBI model. Let us start with the analysis of the peculiarities of the NCBI fluid. 
Looking at the definitions of p nc u and p nc bi (jUj) , <f27l) we can observe that the relation 
between these quantities can be recast in a more significant form: 

Pncbi = Pncbi ~ Pncbi , (30) 

introducing a sort of critical p c ncbi energy density of the model defined as: 

Pncu = /? 2 |r 2 2 (t)(i-rV 2 w)T^[0 ; 0] (i + ^n 2 ^]) 2 } • (si) 

Now, combining with the energy density, the relation (|30"|) can be written as follows: 

Pncbi = (Wncti ~ 1) Pncbi , (32) 

with = y^;- This definition defines the properties of the NCBI fluid. 

One checks easily that the energy density can take on negative values which will lead 
to the so-called ghost solutions. Taking into account the regions of the phase space {<p, <p} 
where the energy density has only positive values we can deduce the conditions on the 
equation of state of the fluid in relation to the value of u>£" 6 *. In fact, depending on the 
sign of this quantity, we shall get different situations as shown in the following table [TJ 



1 1 icfli 

w ncbi 


^Uncbi 


Cosmological behavior 





> 




Inflaton-like 







-1 


Cosmological constant-like 


P 


< 


< -1 


Phantom-like 


]-oo, -/3[U]/3, +oo[ 



Table 1: The cosmological behavior depending on the equation of state of the fluid. 



These results can be easily recovered looking to the definition (p2|). In the last column are 
given the constraints on the scalar field defining the different behaviors. From the defini- 
tion of p c nchi it is easy to observe that the sign of this term and then as a consequence, the 
positivity energy density, of is ruled by the term (1 — j3~ 2 (f> 2 (t)). As a consequence, 
the Inflaton-like rate, the Phantom effect and the cosmological constant-like behavior are 
obtained if respectively <fi is between — (3 and /3, when this quantity is out of this interval 
or if 4> is equal to the NCBI parameter itself. 

Let us now discuss the behavior of this model with respect to the various Energy Con- 
ditions. Again, it will be possible to provide only qualitative and numerical calculations, 
but they will allow us to deduce some interesting properties of the model. We know that 
Energy Conditions can take on several forms. Each of them will define different regions 
in the phase space of the scalar field. 
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The positivity of the energy density parameter is important since it prevents the 
system from possible vacuum instability and rules out the ghost-type solutions. Imposing 
such a condition selects a particular region in the phase space, see Fig. [TJ Moreover, 
the Weak Energy Condition p nc u > and p nc u + p nc u > 0, is strictly related to the 
dynamics of the field since the second relation directly affects the value of w nc u and rules 
out supraluminal solutions. 

The Dominant Energy Condition and the Strong Energy Condition determine some 
areas of the scalar field phase space {0, 0} within which these conditions are satisfied, 
see Fig. El In particular the values of , violating the Strong Energy Conditions lead 
to accelerating solutions. 

The results are plotted in Fig. [T] and Fig. [2] for a particular choice of values of 
parameters (3 and 7, but these plots can be easily generalized to other values of this 
quantities yielding similar or slightly different shapes. The results obtained from the 
imposition of the Strong Energy Condition are interesting, too. In such a case the 
positive value regions are defined around certain values of the scalar field, with large 
areas outside where this condition is violated. This implies that in order to achieve 
accelerated regime in the later stages the scalar field does not necessarily need any fine 
tuning of initial conditions. 



30 
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-20 
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-20 -10 10 20 
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Figure 1: This graph shows the allowed region in the plane {4>(t), <fi(t)} defined by the Weak 
Energy Condition. The plot is obtained in the case of m p i = 10, (3 = 10 7n p ; and 7 ~ m p i. 
The left part shows the energy density curve, the regions inside the curve correpond to negative 
values, which become positive outside. The right panel shows the condition p nc u + PncU > 0, 
which are positive between the two straight lines and negative outside. The = line gives a 
cosmological constant behavior as implied by the Eq. (|31jl. 
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Figure 2: The plot showing the domain in the phase space satisfying the Dominant Energy 
Condition (left) and the Strong Energy Condition (right). The plot corresponds to the choice 
of m p i = 10, j3 = 10 m p i. The left graphic is not influenced by the choice of different values 
for the parameter 7 if (3 is set.h e Strong Energy Condition modifies the contour plots. The 
three different curves in the right part refer to: 7 ~ m p \ (solid), 7 ~ 10 _1 m p ; (dashed) and 
and 7 ~ 10~ 2 m p i (dot-dashed). The inner regions correspond to positive energy values while 
the external ones are negative and provide suitable initial condition for accelerating dynamics. 



4.2 Critical points 

Let us now investigate the dynamical properties of the NCBI scalar field. In terms of 
relevant variables only, the system evolves in a four dimensional space (a(t), a(t), 0, 4>{t)) 
which can be reduced to a three-dimensional set. In fact, in a spatially flat case it is 
possible to introduce the Hubble parameter as an independent variable and to study the 
three-dimensional phase space (H, 0, 0) (here and in the following we neglect the explicit 
time dependence of dynamical variables). Furthermore, the Friedmann equation provides 
a surface in this three dimensional space, implying that the relevant dynamics can be 
followed on this two-dimensional surface. 

However, due to the complexity of the model, it is impossible to get an explicit 
expression of both H and 0. To overcome this problem, a reasonable approach is to solve 
the field equations numerically and to reconstruct the phase space curves for several 
values of parameters. 

The first thing to do while performing a qualitative study is to find the critical points 
of the system. The critical points of the scalar field can be obtained by solving Eq. (f2o|) 
with = and 0(t) = 0, they are obtained the three values = 0, 7v^3/2 or 7\/3. 
On the other side, the critical points of the full system are obtained considering also the 
other field equations. From the cosmological point of view they correspond to de Sitter 
solutions with constant scalar field. The general situation can be summarized as follows: 
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(0, 0) Minkowski 




1 + /TV) 1 ( 1 + ^) ? ) - 1 ) — I <Je Su rer . (33) 




(0,7\/3) Minkowski. 



The first and the third fixed points correspond to trivial Minkowski spaces [321 . In fact 
from the definition of the energy density (|26|) and pressure (|2T|) we see that in these cases 
both do vanish. In other words there is no source for dynamics and the system reduces 
to a Minkowski space-time. 

The second fixed point is a typical de Sitter solution with the equation of state for a 
scalar field of the cosmological constant type. The energy density is different from zero 
and provides a source for the exponential expansion: 

fCH Uer = P 2 1 (1 + /TV) 1/2 (l + ^f) V4 - 1 1 • (34) 

As we will see later, in our model this solution is no more an attractor in the usual 
sense. In fact, there are some values of the parameters for which phase trajectories deviate 
from the de Sitter solution even if the initial conditions are taken very close to this fixed 
point. This result does not agree with the the No-Hair Theorem jHll which holds for 
general cosmological models containing scalar field, as demonstrated in [32] both with a 
minimal and a non-minimal coupling cases. Similar results concerning the dynamics of 
a Born-Infeld scalar field system have been recently obtained in [33J. 

The interpretation of the scalar field fixed points becomes obvious with the Born- 
Infeld relaxed limit analyzed in the previous section. It is easy to check that the two values 
4>{i) = (0, \/37) correspond to the true vacuum configuration of the top hat potential 
in Eq.(J2EJ) while <fi(t) = V3/2'y corresponds to the unstable "false vacuum" state. In 
what follows, we shall sometimes refer to these states as with the usual double- well Higgs 
potential. Let us now analyze the stability of the fixed points. We can rewrite the NCBI 
lagrangian (fT9|) in a new form and rescale the scalar field by substituting 4> — ► \/30: 



C 



ncbi 



1- 1- 



2 



+ 



16 <f) 2 (-7 + 
P 2 



1 + 



4 2 (- 7 + 
P 2 



\2 \ a 



In the vicinity of singular points we observe behaviors similar to those of a scalar field 
minimally coupled to gravity [35J, with energy density given by: 



P = b(-4> 2 + ~m 2 (f> 2 )+p . 



(35) 
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For each fixed point, the parameters m 2 and po should be expressed in terms of the 
parameters (3 and 7 (the parameter b is not relevant for the stability analysis). We just 
saw that for the "Born-Infeld scalar field" there are the three fixed points (|33l) . which in 
terms of the scalar field variables 0,0 are given by (0 = 0, = 0), (0 = 7, = 0) and 
(0=|,0 = O). 

Now, linearizing the equations of motion around these fixed points and evaluating the 
corresponding Jacobian matrices, we infer the exact character of each singular point [36j. 
For the two asymptotic Minkowskian solutions, that is at the points (0 = 0, = 0) and 
(0 = 7, = 0), we find that 

m 2 = 4 7 2 po = . (36) 

which corresponds to a stable point with eigenvalues in the (0, 0) plane: A± = ±^7, 
both complex. In the case of the point (0 = |, = 0) we find 

m . = _v-^- P0 = (i + i;) , v /r^Li. (37) 

Hence, because m 2 is negative, it corresponds to an unstable point, and the eigenvalues 
in the plane (0, u) are: 

3# , / 4^ , , 3H 0/ / 4m^ N , . 

X + = ^f-9Hl " 1)>0 A - = "T^ " V 1 - 9sf < ° < 38 > 

with Hq = J yPo- In terms of (3 and 7, the eigenvalues take on the following form: 



3 




\ 



-1 + 



\ 



32 7 2 (2/5 2 + 7 4 ) 



3 (4/3 2 + 7 4 ) 



-2+H + 



4 + ? 



(39) 
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A_ 



-3 
2V$\ 



( 



1+ 



V 



\ 



32 7 2 (2/3 2 + 7 4 ) 



3 (4/5 2 + 7 4 ) 



1 + ? 



4 + S 



(40) 



/ 



The most interesting feature for the stability study is the sign of A + and A_ which 
is not easy to determine. However, after the expansion in powers of (3 for A± these 
expressions are considerably simplified: 



A, 



^7- 



7 4 1 



-V2' 



7 -4\/2^ + ° ( ^ ; 



(41) 



4.3 Numerical Analysis 

Our dynamical system is characterized by three parameters (3, 7, mpi, which in the 
case of flat space geometry determine it completely. However it is obvious that only two 
of them are independent. To perform the numerical analysis we have considered /3, 7 
in terms of the normalized Planck mass. In order to explore a dimensionally uniform 
parameter space we consider 1/(3, 1/7 2 , \ jm 2 pl 3 . 

The parameters (3 and 7 cannot take on arbitrary values if we have to satisfy the 
following physical requirements: 

i) the characteristic Born-Infeld parameter cannot be higher than the squared Planck 
length, i.e. (3 < m 2 Pl ; 

%%) the mass parameter also can not be greater than the Planck mass, i.e. 7 2 < m Pl . 

Taking into account these constraints we have explored the space of the parameters 
considering a domain ranging from fractions of the Planck constant to its entire value. 
To perform a numerical analysis we have fixed the value ofmp; at 10, so that the essential 
features of the phase space behavior can be clearly displayed. To give a picture of the 
results we have considered (3 G [10~ 2 m Pl ,m Pl ] (namely with mpi = 10, 1/(3 G [10 -2 ,1]) 
and 7 G {lO- 2 m Pl , m Pl ] (I/7 2 G [10 -2 , 100]). 

The definition of the lagrangian and the field equations impose some natural con- 
straints in the phase space {0 ,<p}. The lagrangian (fTTl contains an even root term, this 
means that the term under the root must be positive. Indeed it is as it is a sum of two 
squared terms. Let us look at the field equations. The Friedmannn equation (|23|) . in the 
spatially flat case we are considering, implies the positivity of energy density p nc u > 0. 

3 We remember that the Planck mass has the dimension of [cm] -1 , while the Born-Infeld parameter 
has the dimension of [cm] -2 and the mass parameter 7 has the dimension of [cm] -1 . 
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This gives the first constraint on the phase space of the scalar field as we already saw 
in FigHJ Let us also note that there are no divergences in the definition of the energy 
density and that it remains always finite. 

Considering now the pressure equation we observe that the only constraint is 
again the positivity of the function under root, but it leads to the same constraints as 
those imposed by the lagrangian itself. 

The analysis of the scalar field equation (|25|) is more complex. Here we observe that if 
the function multiplying the (j)(t) term is zero than the term (j)(t) can become divergent. 
This implies, as in the case without gravity (Hj, that there exists a characteristic curve in 
the phase space where the solutions for the scalar field become divergent in a finite time. 
This problem can be cured if all other terms in the equation simultaneously vanish. This 
allows (f)(t) to remain finite and provides some crossing points in the phase space on the 
separatrix were the system can go beyond this curve. The singular curve 'separatrix) is 
given by the following relation: 

(1 + V[0(t)]) { (1 - /r 2 2 (t)) 2 + yn 2 [0(t)] (1 - 3/r 2 2 (t)) } = o . (42) 

Now, comparing with the pure Higgs field case studied in [H] there is an obvious and 
important difference. The presence of gravity, which induces the new constraint p nc u > 
implies the absence of crossing points over the divergence curve. 

With the allowed regions thus defined, the following step is to plot the characteristic 
curves generated by the scalar field equation. The cosmological evolution of the model 
and the scale factor behavior are driven by the scalar field dynamics. 

At this point we are ready to perform a numerical study of the system. In the follow- 
ing we will present the results for three combinations of parameters expressed in terms 
of the Planck mass: 

Case: (3 = m 2 Pl . 

The system displays similar phase trajectories for all values of mass parameter (we plot 
the case 7 ~ mpi and 7 ~ 10~ 1 mpi). The scalar field shows a spiralling behavior ending 
up in one of the two stable fixed points. 

The time evolution on phase trajectories is from left to right. For smaller values of 
7 the allowed region around the fixed points becomes narrower. This is understandable 
because the mass parameter in the standard 4 (t) theory controls the width of the poten- 
tial. One obvious difference with the case analyzed in [Bj is the spiraling behavior. This 
effect is due to the presence of a friction term in the scalar field equation coming from the 
gravitational interaction. The pure Born-Infeld regime is recovered when one increases 
the values of the Planck mass (equivalent with a decrease of gravitational coupling) and 
fixes the other parameters. On the other side if one increases the (3 parameter the effect 
is to widen the allowed domain along the <j)(t) direction. This transforms the system into 
the standard Higgs case in presence of gravity obtained in the relaxed limit (fl~8|) . 
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Figure 3: Phase space of the scalar field in the case (3 = m pi with 7 ~ mpi and 7 ~ 10 _1 mp;. 
The dashed curve is the constraint obtained by the Weak Energy Condition (as yet seen in 
Fig^), while the lighter one is the characteristic curve due to the scalar field equation. We show 
only some characteristic curves since for other initial conditions the shape is analogous. 

A particular result visible in Fig.(|BJ is that the physically significant dynamics is 
constrained to a very small region around the fixed points. All the trajectories with 
different initial conditions intersect the singular curve in a finite time. This entails a 
singular behavior for the scale factor thus excluding such solutions. It is easy to find 
the cases in which it is possible to get inflationary behavior. In general (see TabEJ) 
non-linearities combined with a small enough mass parameter make the scalar field roll 
down too quickly, never attaining a slow rolling regime. This implies that the scale factor 
follows a power law regime since the beginning. 

This conclusion results from examining the slow rolling parameter e = — |29| . It 
is well known that this parameter should be significatively smaller than 1 throughout all 
the inflationary expansion. 

In the cases 7 ~ rn Ph 7 ~ 10~ 1 mp; and 7 ~ 10 _1 mp; with a slow rolling-type initial 
condition from the top of the unstable fixed point, it is possible to get inflationary 
evolution with the right power law phase transition at a given stage. 
Case: (3 = m 2 Pl /10. 

We examine the plot with 7 of the same order of the Planck mass. This plot is quite 
interesting since it is clear that the two stable fixed points define two separate regions 
between which there are no trajectories. The other cases with smaller values of the mass 
parameter recall quite similar shapes to the previous case. 

Around the two fixed points there are stable spiralling solutions. The trajectories starting 
from the unstable fixed point hit the singular curve in a finite time. They describe closed 
lines broken by the characteristic curve. If the initial velocity is small enough, the scalar 
field remains in the unstable critical point for a long time and then it joins again the 
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Figure 4: Phase space of the scalar field in the case (3 =,m Pl /10 with 7 ~ mpi. The dashed 
curve is the constraint imposed by the Friedmann equation, while the thinner one is due to 
the scalar field equation. We limit the number of solutions plotted to make the results clearly 
visible. 

characteristic curve in finite time since the remaining part of the phase space is forbidden. 
This provides an infinite inflationary solution of cosmological constant type. It is obvious 
that in this case the de Sitter solution is no more an attractor. The initial conditions 
taken close enough to the unstable fixed point instead of providing inflationary de Sitter 
like solutions give power law trajectories which never collapse into the stable vacuum 
state intersecting the = 00 curve instead. 

In other cases the behavior is similar to the previous case with (3 = m Pl . Again 
no inflation is observed, and the slow rolling parameter e is always greater than one. 
On the other side the phase space obtained with 7 = 10~ 1 mpi provides the physical 
condition ensuring inflation and it remains similar to the one shown in Fig.Q. All 
possible cosmological behaviors of the model are given in Tab.(|2|), also for the values of 
parameters not considered in the plots. 

We have tested also the amount of inflation provided by the NCBI model. Via nu- 
merical evaluation of it was possible to calculate the winding number of trajectoriy 
by means of the relation 



We see that the NCBI scalar model can generate the right amount of inflation (more than 
60), if proper slow rolling conditions are imposed. The expansion rate of this model is 
lower than the one obtained with the corresponding Born-Infeld relaxed lagrangian (|29|) . 
Case: (3 = 10" 2 m^. 

This case represents the strong Born-Infeld limit in relation to the small value of the 
Born-Infeld parameter. The shape of the phase space is similar to the ones proposed 
before, with the big feature that the allowed region is strongly constrained. From the 
cosmological point of view the behavior is quite the same of the case with f3 = 10 _1 mp ; . 




(43) 



21 



mpi = 10 


(3 -> 10m^ 




io- 1 ™^ 




l(T 3 m^ 




10- 5 m 2 Pl 


7 -> m P i/V3 


— 


+ 


+ 


+ 


+ 


— 


— 


T ^= 




+ 


+ 


+ 


+ 


+ 




m Pl lO- 2 /y/3 




+ 












m Pl lO- 3 /V3 
















m P ilO- 4 /V3 
















m P ilO- 5 /V3 
















mpilO- 6 /V3 

















Table 2: A summary of the capability of providing inflation for the NCBI cosmological model. 
The combination of parameters allowing inflation are indicated with a sign +, the wrong ones 
whit a — . It is obvious that this model can provide inflation only if the mass of the scalar field 
is close the Planck mass. With small values of (3 we get the strong BI regime, it appears that 
in this case only some particular configuration of parameters allow inflation. 

In relation to the variation of 7 is possible to obtain inflation only for some condition of 
the mass parameter. In FigElwe show two cases, again the solution 7 = 10~ 1 mpi shows 
a particular behavior as in precedence, the following cosmological properties are similar 
to the previous case. Another interesting phenomena appears even in the second plot 
when the mass-shell of the scalar field is close to the Planck scale. In such a case the 
wrapping of the phase lines evolves moving around both the two static fixed points while 
the stationary configuration is achieved ending in one of this points in relation to the 
initial condition without any particular rule, see Fig|3 This happens for initial conditions 
which are chosen for all the three fixed points with several values of the velocity 0. 
Interesting plots are shown in Figlll In this case a configuration of the parameters 
is chosen in such a way that the gravitational coupling happens to be relaxed while 
the BI coupling is strong. The net effect is that the spiralling behavior around stable 
fixed points in the left plot of Figgis more similar to the one proposed in Closed 
trajectories shown in absence of gravity are no more possible due to the friction term 
and are transformed into spirals. No significant differences are found with the variation 
of the mass parameter. 

The simultaneous evolution of the cosmological scale factor a(t) (through its logarith- 
mic derivative with respect to In t) and of the scalar field is displayed in Fig. [HJ 
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7 = W~ 2 mpi. In the second plot the region in the middle of the curve trajectories (small values 
of (f), left in white) has been plotted in detail in figure E3 It is obvious a behavior similar to 
the one found in jS] affected by the spiraling shape of the trajectories due to the gravitational 
interaction. 




* 



Figure 6: The detail of the previous plot, it is showed how the wrapping of the solutions 
in proximity of the critical points changes in relation to the initial conditions, left side (ft = 
0, (j) = 0.16 and right § = 0, <fi = 0.18. This behavior is alternated, for cj) = 0.36 the 
solutions converges to <f> = while for <j) = 0.6 again versus eft = 0.33. 
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Figure 7: The phase portrait in the case of low gravitational coupling. A strong similarity is 
found with the results of [S]. The gravitational interactions break the steady state configuration 
of solutions; the thickness of the line is due to the very slow spiralling of the trajectory. 
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Figure 8: The simultaneous evolution of the Hubble parameter and of the scalar BI field. The 
last graph is the phase space plot of (f>, <j>; the trajectory is spiralling towards the singular point. 

5 Discussion and conclusions 

In this paper we propose the cosmological analysis of a NCBI model obtained from the 
matrix realization of a non-commutative geometry. This study represents the gravity 
coupled generalization of a model proposed in [6j. To determine the cosmological proper- 
ties of the model and its physical relevance in this sense we have investigated numerically 
the NCBI scalar field dynamics constructing its phase space characteristic curves. The 
cosmological results, in relation to several values of the parameters, are summarized in 
Tab.©. 

It is possible to observe that there are different behaviors, but in any case inflation 
can be recovered when the mass of the scalar field is close to the Planck mass. A relevant 
difference with the 4 (t) limit, is that the obtained amount of inflation is smaller, although 
enough to solve the several shortcomings of the standard cosmological models. 

The scalar field shows a phase space which is similar to the standard case, without 
a BI cut-off, if (3 is big in comparison to the other parameters. In the other cases the 
dynamics develops inside a constrained region around the critical points. However, the 
bounded regions distinctive of the Born-Infeld free theory pa] are lost because of the 
presence of the cosmological friction term. A shape similar to the pure BI case can be 
observed in Fig. but it is obvious a slow spiralling rate of the trajectories as a slowly 
decaying in the stable vacuum minima. 

The most important difference with the standard case is that in the BI realm all 
the dynamics is always constrained inside a well defined region. This is a persisting 
feature which can be found in all phase space plots. Since for a strong BI coupling the 
physically significant region is very reduced both in the velocity and in the scalar field 
values, this means that such a framework possesses an intrinsic condition that the scalar 
field should move starting only near the critical points. This is more obvious when the 
Born-Infeld parameter is small. All other initial conditions will provide trajectories which 
will intersect the singular line in a finite time leading to physically meaningless solutions. 
In some sense this approach contains in itself a particular choice of initial conditions 
without any fine tuning procedure. 

In order to decide whether the scalar NCBI model is a viable approach of inflation, still 
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a large amount of work is needed, above all with respect to the observational predictions 
of the model. It should be interesting to verify whether such a model can represent a 
good candidate for the quintessential inflation (3Zj. It is interesting to ask if the high 
non-linearity of the NCBI approach can preserve safe contributes even in the late time 
epoch able to generate a driving force for the new accelerating phase of cosmic expansion. 
These questions remain beyond the scope of the present paper. 
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